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A precision measurement of the p{e,e'p)7r^ reaction near threshold 
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New results are reported from a measurement of tt® electroproduction near threshold using the 
p(e, e'p)7r® reaction. The experiment was designed to determine precisely the energy dependence 
of s— and p—wave electromagnetic multipoles as a stringent test of the predictions of Chiral Per¬ 
turbation Theory (ChPT). The data were taken with an electron beam energy of 1192 MeV using 
a two-spectrometer setup in Hall A at Jefferson Lab. For the first time, complete coverage of the 
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and angles in the p7r° center-of-mass was obtained for invariant energies above threshold 
from 0.5 MeV up to 15 MeV. The 4-momentum transfer coverage ranges from 0.05 to 0.155 
(GeV/c)^ in fine steps. A simple phenomenological analysis of our data shows strong disagreement 
with p—wave predictions from ChPT for > 0.07 (GeV/c)^, while the s—wave predictions are in 
reasonable agreement. 

PACS numbers: 25.30.Rw, 13.60.Le, 12.39.Fe 


Neutral pion production from the proton vanishes in 
the chiral limit of zero quark masses and pion momenta 
Ptt —f 0. As a result, the reaction at threshold is partic¬ 
ularly sensitive to non-perturbative mechanisms within 
QCD which break chiral symmetry. It is also experimen¬ 
tally the most challenging to study. Pion photo- and elec¬ 
troproduction experiments are now producing data of un¬ 
precedented precision to test Chiral Perturbation Thec^ 
(ChPT), the low-energy effective field theory of QCD [l|. 
ChPT treats the spontaneous and explicit chiral sym¬ 
metry breaking in terms of a perturbative expansion in 
small momenta and quark masses, and makes predictions 
for the s— and p—wave multipoles for the yiV —>■ 7riV re¬ 
action in the near-threshold region. Within ChPT, the 
internal structure of the pion and nucleon is systemat¬ 
ically parameterized by Low Energy Constants (LEC), 
while the long-range external "kN dynamics are fixed by 
the underlying chiral symmetry. Once the LECs are de¬ 
termined by experiment near threshold, the convergence 
of the chiral expansion can be tested by comparing pre¬ 
dictions with data taken at energies above threshold. 


Recently 7r° photoproduction cross-section and po¬ 
larized photon beam-asymmetry (E) data from the 
MAMI A2/CB-TAPS experiment |3| were used to test 
two versions of ChPT. The relativistic ChPT calcula¬ 
tion (RChPT/yMAID) has been carried out to 

O(p^), while the non-relativistic Heavy Baryon ChPT 
calculation (HBChPT) is of 0{pj^) for photoproduction 
(BKMOl) [6| but only of O(p^) for p—waves in electro¬ 
production (BKM96) 0. Both the BKMOl and RChPT 
calculations, after fits of LECs to the data, were com¬ 
patible with the experimental multipoles Ao+j Ei^, 
and Ml- within an incident photon energy range of 7 to 
25 MeV above threshold [3, [8|. 


The pion electroproduction reaction 7 *p p7r° allows 
a more stringent test of ChPT, since the four-momentum 
transfer and invariant energy W can be varied in¬ 
dependently. Chiral ttN dynamics naturally involve the 
mass scale I while the LECs fitted in photoproduc¬ 
tion encapsulate higher order processes, involving pos¬ 
sibly TV A or p,uj degrees of freedom. The depen¬ 
dence near threshold may reveal the onset of these short- 
ranged mechanisms. Until now, only limited kinematic 
coverage from 7 *p —> p7r*^ threshold experiments is avail¬ 
able |9l4l2l|. Several older MAMI experiments showed 
a dependence of the total cross section near thresh¬ 
old incompatible with HBChPT 11, 1^, although a new 
MAMI re-measurement has superceded those data 0. 


The JLAB/Hall A experiment reported here provides the 
most extensive (Q^, W) coverage of 7r° electroproduction 
to date for testing theories of chiral dynamics substan¬ 
tially above threshold. 

Under the one-photon-exchange approximation, the 
p{e, e'p)TT^ cross section factorizes as follows: 

TV 

dQ'^dWdn* ^ d^ ^ ’ 

where r„ is the virtual photon flux and the Jacobian 
J = d(Q^,W)/d(Ee>,cos6e') relates the differential vol¬ 
ume element of data binned in dQ^dW to the scattered 
electron kinematics dEe'd cos 6*e'. The pir^ center-of- 
mass (C.M.) differential cross section, dcr/dH*, depends 
on the transverse e and longitudinal polarization of the 
virtual photon through the response functions Rt,Rl 
and their interference terms Rlt and Rtt'- 

7 * 

^ {Rt + €lRl + e Rtt cos 2(/)* 

-f y/2eLie + l)i?LTCOS(/*). (2) 

The response functions depend implicitly on W and 
0*, the 7r° C.M. angle, while the angle (j)% defines the 
rotation of the pir^ plane with respect to the electron 
scattering plane (e,e'). Other definitions are = 
((5^/|A:*p)e, r„ = aEe>k*W/2Tr‘^Ef.mpQ^{l — e) and 
J = TTW/EgEgimp. Finally |fc*| andp* are the C.M. mo¬ 
menta of the virtual photon and pion respectively, while 
k* = {W^ — mp)/2W is the real photon equivalent energy. 

The p(e, e'p)Tr^ experiment was performed in Hall A 
at Jefferson Lab using the Left High Resolution Spec¬ 
trometer (LHRS) 0 to detect the scattered electron 
and the BigBite Spectrometer m to detect the coinci¬ 
dent proton. The CEBAF beam was energy-locked to 
1192 MeV and delivered to a 6-cm long, 2.54-cm wide 
cylindrical liquid hydrogen (LH 2 ) target. Beam cur¬ 
rents below 5 pA were used to limit the singles rates in 
both spectrometers. Four angular settings for the LHRS 
{9e' = 12.5°, 14.5°, 16.5° and 20.5°) covered a nearly con¬ 
tinuous range of 0.05 — 0.155 (GeV/c)^ using a 4.4 
msr acceptance cut. The LHRS momentum acceptance 
was centered on the p7r° threshold and covered the range 
—3% < Sp/p < -1-5%. 

Three angular settings of the BigBite were used 
(0p=43.5°,48° and 54°) which provided full coverage 
(Fig. [T|) of the proton cone up to an invariant energy 
above threshold of A1U=15 MeV (at the largest Q^). 
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FIG. 1. Left: Overlap between three BigBite Spectrom¬ 
eter proton laboratory angle settings (colored boxes) and 
p7r° center-of-mass bins at Q^=0.135 (GeV/c)^ and AVF=9.5 
MeV. Radial and concentric lines separate bins of and 6’^, 
respectively. Only 5 out of 9 9’^ bins are shown. The blue 
line shows 0^=180. Right: Radial and concentric lines sepa¬ 
rate bins of 9’^ and AW, respectively, projected onto proton 
lab momentum pp and 9p. Bins to the (left,right) of the blue 
line correspond to (0^=180°, (())(=0°). The innermost circle 
represents AW = 0.5 MeV. 


The BigBite momentum acceptance covered the range 
(0.25< Pp <0.5 GeV/c), limited by the target energy loss 
at low momentum and the thresholds on the E—AE scin¬ 
tillator counters at high momentum. The low momentum 
cutoff was achieved using a thin (25 pm) Ti exit window 
in the target scattering chamber and a helium bag for 
transport up to and between the BigBite drift chambers. 
Absolute normalization, energy and angle calibrations in 
both spectrometers were checked at each kinematic set¬ 
ting using elastic scattering runs with LH 2 and thin solid 
targets. 

Scintillator hodoscopes provided the primary triggers 
for both spectrometers. A gas threshold Cerenkov de¬ 
tector in the LHRS provided electron identification with 
99% efficiency. Signals from either E or AE scintillator 
planes at the rear of BigBite were used in the coincidence 
trigger, while signal thresholds in both the hodoscopes 
and multi-wire drift chambers were set to suppress mini¬ 
mum ionizing tracks from pions. Final proton identifica¬ 
tion was made using E—AE cuts on the highly segmented 
scintillators. The path-length corrected coincidence time 
distribution between the LHRS and BigBite is shown in 
Fig. A 10 ns wide cut centered on the peak was used 
to select true coincidences, while a 30 ns cut (excluding 
the peak) selected random coincidences for subtraction. 
Selection of the p7r° final state required calculation of 
the missing-mass M after reconstruction of the detected 
particle’s 3-momenta: 

= {E + mp- Ae/ - £'p)2 - (Pe - Pe' - PpY (3) 
The experimental missing-mass distribution is also shown 




COINCIDENCE TIME (NS) MISSING MASS - Mn° (MEV) 


FIG. 2. Left: Coincidence timing between the LHRS and 
BigBite. Events belonging to the true coincidence peak were 
selected using cuts indicated by the vertical lines, while ran¬ 
dom coincidences were selected from the region highlighted in 
red. Right: Missing mass distribution at Q'^=0.15 (GeV/c)'^ 
for the invariant mass range 0 < AW < 10 MeV. Background 
events from random coincidences (red) and target cell win¬ 
dows (blue) were subtracted from the raw distribution, leav¬ 
ing the 7r° missing mass peak shown in gray. 


in Fig. [5] before and after subtraction of both random co¬ 
incidences and target-window contributions. The latter 
background was estimated using cuts on AW below the 
7r° threshold. 

Before binning the data, both incident and scattered 
electron energies were corrected for ionization losses in 
the LH 2 and target windows, using the calculated en¬ 
trance and exit paths with respect to the measured tar¬ 
get interaction vertex. Proton transport energy losses 
through the target, Ti window and BigBite were also 
corrected for each event. Acceptance corrections were 
derived from a Monte-Carlo simulation of both spectrom¬ 
eters, using the Dubna-Mainz-Taipei (DMT) model [l^ 
as a physics event generator. Special care was taken to 
incorporate into the simulation radiative correction and 
straggling losses, a fine-mesh magnetic field map for the 
BigBite, and the measured energy and angular resolution 
and energy calibration determined from elastic scattering 
runs, in order to properly account for their systematic ef¬ 
fects near threshold. The dominant sources of systematic 
uncertainty are target window background subtraction, 
accidental coincidence corrections and LHRS central mo¬ 
mentum calibration, which combined contribute to the 
overall normalization error of 20% near threshold at low 
decreasing to 7% for data above threshold at higher 
Q^. 

Events were accumulated using (12,30,18,9) bins for 
(Q^, AW, 0*) respectively, with a cut of ±10 MeV on 
the missing-mass peak. The AW bin width was 1 MeV 
and the LHRS acceptance extended up to A1F=30 MeV, 
although with reduced C.M. coverage. The average 
bin width was 0.01 (GeV/c)^. Figure [3] shows typical dif¬ 
ferential cross sections for each (/)* and d* bin obtained at 
(5^=0.135 (GeV/c)^ and AW = 9.5 MeV. The curve la- 
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FIG. 3. Differential cross sections for p(e, e'p)Tv^ from this 
experiment at = 0.135 (GeV/c)^ and AW = 9.5 MeV 
binned in pn^ center-of-mass angles and and cos 6*. See 
text for description of curves. Units are phjsr. Errors are 
statistical only. 
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FIG. 4. Total cross section for p(e, e'p)Tv° as a function of 
for different bins in AW (invariant mass above threshold) for 
(□) this experiment and (A) MAMI [T^. Units of AIU are 
MeV. Errors are statistical only. 


beled BKM96 is the HBChPT prediction from Bernard et 
al. 0, which uses LECs fitted to older photoproduction 
data from MAMI and electroproduction data at Q^=0.1 
(GeV/c)^ from MAMI and NIKHEF. The other curve is 
an empirical fit to the data which we use to obtain the 
total cross section atot- The empirical fit uses the form in 
Eq.(2) and expands the response functions with Legendre 
polynomials Pi{x), where x = cos0*, 

Rt + ^lEl = Al+^ + Al+^ Pi{x) + Al+^ P2{x) (4) 
Rtt = A^ (1 - x^) (5) 

i?LT = + Af'^Pi(x)) (1 - (6) 

The total cross section atot is given by 47rf5- Aq~^^ . 

The dependence of atot is shown in Fig. H] for dif¬ 
ferent AW bins starting 0.5 MeV above threshold. Two 
ChPT calculations are shown (BKM96 [ 3 ], x^AAID [^), 
along with the SAID08 solution 119 and phenomeno¬ 
logical models (DMT [^, MAID ^3|) which have been 
fitted to the world data on pion photo- and electropro¬ 
duction. Compared to the linear dependence of the 
HBChPT/BKM96 curve, our atot measurement shows a 
bending over at higher similar to the phenomenolog¬ 
ical models and the RChPT/xMAID theory. At lower 
Q^, both ChPT calculations are consistent with our data 
over the entire AW range shown here. Note that two 
of the RChPT LECs were fitted to a new MAMI re¬ 
measurement [l^ (triangles in Fig. 01) of earlier > 0 
experiments, while the remaining LECs were fitted to the 
= 0 A2/CB-TAPS data 0. 

Near threshold, the s— and p—wave decomposition of 
atot can be obtained by fitting the p* dependence of 
Aq according to 

A'^+^ = ao + b\pl\^. ( 7 ) 




FIG. 5. The dependence of ao (left) and b (right) from the 
fits of Eq. 7 to the Legendre coefficient . The theory 

curves are calculated for the beam energy of our experiment 
(1192 MeV). For the curve labeled REFIT the BKM96 LEG 
bp has been lowered from 13 to 9.3 (GeV)~® (see text and 
Fig. 6). Errors are statistical only. 


The b coefficient parameterizes the contribution of 
p—wave multipoles arising from their intrinsic p* de¬ 
pendence near threshold, while ag fits the combination 
|Eo+P -I- eilLo+P of s—wave multipoles extrapolated to 
threshold. The Lo+ multipole dominates ag over our 
range due to a large ep factor. The extraction of ag and 
b from fitting our data up to AW = 9.5 MeV is shown 
in Fig. [5l along with fits to the newest MAMI data [l^ 
up to AW = 3.5 MeV (the limit of their measurement) 
and previous results from NIKHEF There is good 

agreement of both ag and b with the chiral model pre¬ 
dictions for our lowest points. For higher the 
HBChPT curve describes ag better than RChPT. How- 
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FIG. 6. The AW dependence of Legendre coefficients from 
fits to our data at = 0.135 (GeV/c)^. Note that the 
AVF=9.5 MeV bin corresponds to the Legendre fit shown in 
Fig. [3] Errors are statistical only. 


ever the strong disagreement of our b coefficient with both 
chiral curves for > 0.07 GeV^ suggests at least one 
of the p—wave multipoles is described incorrectly in the 
calculations. The dependence of b from fitting the 
MAMI data is qualitatively similar, although with larger 
errors, due to the smaller AW range of their data. 

Further insight can be obtained from the AW de¬ 
pendence of the Legendre coefficients in the > 
0.07 (GeV/c)^ region. This is shown in Fig. | 6 ] at 
Q^=0.135 (GeV/c)^. While all models are in good agree¬ 
ment with our data near threshold, the theory curves 
for and Aq^ show large variations above 

AVF=3 MeV. These coefficients are particularly sensitive 
to the p—wave multipole combinations P 3 = 2Mi+-|-Mi_ 
and P 2 = 3Pi+ — Mi_|_ -|- Mi_, while Aq ^ is also sensitive 
to the combination APfg = (Pf “ Our fit result 

for is close to zero over the ALF range of our data, 
which implies P| ~ Pf or ss —2 (neglecting 

the weak electric quadrupole Pi+). 

Only the DMT model predicts k, 0 for AW < 
15 MeV, largely due to their calculation of the Mi_ mul¬ 
tipole [l9| , the value of which is substantially larger than 
predicted by ChPT. A similar result was obtained from 
disperson relations [i^. In the BKM96 theory, which 
uses a 0{p^) p—wave expansion, it is not possible to sep¬ 
arately adjust A/i+ and A/i_, since only P 3 is controlled 
by a single LEG bp. By reducing 6 p in the calculation 
from 13.0 to 9.3 (GeV)“^, we can improve agreement 
with both A^^^ and Aq ^ as shown in Fig. | 6 ] by the 
curve labeled REFIT. However this adjustment worsens 
the agreement with p—waves at lower as indicated 


by the REFIT b curve in Fig. [5] Moreover, a different 
adjustment oi bp is required to match our measurement 
of Aj+^. 

The 0{p^) RChPT calculation i predicts a nearly 
identical dependence for the b curve in Fig. [S] as the 
0{p^) HBChPT theory. At leading-order and next-to- 
leading order, P 3 is controlled by a single 0{p^) LEG 
dg, similarly to HBGhPT [^. However dg is highly con¬ 
strained by the < 3^=0 photoproduction fits, and there 
is almost no room for adjustment. Other 0(p'^) LEGs, 
which explicitly control dependent terms, either do 
not appreciably affect the p—wave multipoles, or effect 
the same response as bp. 

Despite the very different LEG composition of 
HBChPT and RChPT, it appears neither calculation can 
be adjusted to agree with the trend of our p—wave 
data. Purthermore this discrepancy occurs well within 
the AW range wheryrhotoproduction p—waves are well 
described at 0 {p ^) Q . Our data therefore suggest that 
higher powers of are needed in the ChPT formalism, 
while the onset of disagreement {Q^ >0.07) implies a 
t—channel energy scale above the pion mass. Similar dis¬ 
crepancies in ChPT calculations of nucleon form factors 
were removed by including vector mesons as dynamical 
degrees of freedom [2l|. Our data could provide strong 
constraints to analogous extensions of pion electropro¬ 
duction calculations. 

In summary, a JLAB/Hall A experiment has measured 
for the first time both the and extended AW depen¬ 
dence of the threshold p(e, e'p) 7 r° reaction with full C.M. 
coverage and fine binning. Our phenomenological fit of 
the data shows reasonable agreement with two leading 
ChPT theories for s—waves, while chiral predictions of 
p—wave contributions strongly diverge from our data for 
> 0.07 (GeV/c)^. We use a Legendre decomposition 
of our cross sections to show there is insufficient flexibil¬ 
ity in the low energy constants available for p—waves to 
account for the discrepancy. 
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